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Abstract 

Radiative returns to the Z resonance (Z7 events) are used to determine the LEP2 centre- 
of-mass energy from the data collected with the ALEPH detector in 1997. The average 
centre-of-mass energy is measured to be: 

E CM = 182.50 ± 0.19 (stat.) ± 0.08 (syst.) GeV 

in good agreement with the precise determination by the LEP energy working group of 
182.652±0.050 GeV. If applied to the measurement of the W mass, its precision translates 
into a systematic error on My/ which is smaller than the statistical error achieved from 
the corresponding dataset. 
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1 Introduction 



One of the main goals of LEP2 is the direct measurement of the W mass with an accuracy 
better than 30 MeV. This requires a very precise determination of the LEP centre- 
of-mass energy, in order to minimize its contribution to the W mass uncertainty. At 
LEP1 this accuracy was achieved using the Resonant Depolarisation method for energy 
measurement [jlj. At LEP2 this is no longer possible since polarization has not been 
achieved above a beam energy of 60 GeV. Therefore, the precise Resonant Depolarisation 
calibration is performed only in the 40-55 GeV energy range, and extrapolation to higher 
energies is necessary. The precision achieved on the 183 GeV centre-of-mass energy range 
is about ±50 MeV 0. For future data taking there is hope to improve this precision 
using the recently installed magnetic spectrometer || . 

An alternative determination is considered in this paper. It is based on an analysis of 
e + e~ — > 7// events, where the Initial State Radiation (ISR) photons are mostly collinear 
to the beam pipe and remain undetected. In these events, the invariant mass, a/s 7 , of 
the final state // system peaks naturally at the Z mass. The variable x — 1 — s'/s can 
be reconstructed from the final state particle angles which are well measured quantities. 
Since the Z mass is very well known from LEP1, a fit to the distribution of x can be used 
to determine the center-of-mass energy y/s. 

This letter presents an application of this method using the qq final state in data 
recorded by ALEPH in 1997. LEP ran at nominal centre-of-mass energies of 181, 182, 
183 and 184 GeV, and an integrated luminosity of 56.812 pb" 1 was collected. 

2 The ALEPH detector 

A detailed description of the ALEPH detector can be found in Ref. || and of its 
performance in Ref. Hj. Charged particles are detected in the central part of the detector. 
From the beam crossing point outwards, a silicon vertex detector, a cylindrical drift 
chamber, and a large time projection chamber (TPC) measure up to 31 coordinates 
along the charged particle trajectories. A 1.5 T axial magnetic field is provided by a 
superconducting solenoidal coil. A resolution of Spx/pr = 6 x 10~ 4 j>r © 0.005 (p? in 
GeV/ c) is achieved. Hereafter, charged particle tracks reconstructed from at least four 
hits in the TPC, having a polar angle to the beam axis satisfying |cos#| < 0.95, and 
originating from within a cylinder of 2 cm radius and 20 cm length, centred on the 
nominal interaction point and parallel to the beam axis, are called good tracks. 

Electrons and photons are identified in the electromagnetic calorimeter by their 
shower profile. The calorimeter, a lead/ wire-plane sampling device with fine readout 
segmentation and total thickness of 22 radiation lengths at normal incidence, provides an 
energy resolution AE/E of 0.180/v 7 ^ + 0.009 (E in GeV). 

Muons are identified by their penetration pattern in the hadron calorimeter, a 
1.2 m thick iron yoke instrumented with 23 layers of streamer tubes, together with two 
surrounding layers of muon chambers. In association with the electromagnetic calorimeter, 
the hadron calorimeter also provides a measurement of the energy of charged and neutral 
hadrons with a relative resolution of 0.85/a/E (E in GeV). 
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The total visible energy and momentum, as well as the missing energy, are evaluated by 
an energy flow reconstruction algorithm || which combines all of the above measurements, 
complemented at low polar angles by the energy detected in the luminosity calorimeters. 
Jets are built from charged and neutral objects reconstructed by the energy flow algorithm. 
The typical jet angular resolution is 30 mrad. The jet energy resolution is approximately 
<je = (0.6y/E + 0.6)GeV ■ (1 + cos 2 ^), where E (in GeV) and 9 are the jet energy and 
polar angle. The jet energy and angular resolution as well as calibrations were obtained 
from extensive studies of Z — > qq events both in data and Monte Carlo. Discrepancies 
between data and simulation were used when evaluating systematic errors. 



3 Monte Carlo samples 

A sample of 200,000 qq events was generated using KORALZ v4.2 at a nominal energy 
of 182.675 GeV and fully simulated in the detector. Monte Carlo samples at 183 GeV 
with integrated luminosities at least 20 times larger than recorded were simulated for all 
background reactions. PYTHIA v5.7 [§] was used to generate ZZ and Zee events, and 
KORALW vl.21 |J to produce four fermion events with WW topologies. Two-photon 
(77) reactions into hadrons were simulated with PHOT02 [11|| and PHOJET ||12|| . The 



method was calibrated with fully simulated qq samples of 20,000 events, generated with 
KORALZ, at 181, 182, 184 and 185 GeV. 



4 Event selection and reconstruction algorithm 

At t/s = 183 GeV the main backgrounds to the process e + e~ — > Z7 — > qcp( are 
e + e~ — > WW, e + e~ — > ZZ, e + e~ — > 77 and e + e~ — > Zee. To discriminate between 
signal and background events, the following cuts are applied. 

• qq candidates are required to have at least seven good tracks which total energy 
exceeds 10% of the nominal centre-of-mass energy. 



The energy flow objects are clustered into jets using the JADE algorithm |T3| with 



a y C ut value of 0.008. Events containing jets with an electromagnetic content of at 
least 90% and an energy larger than 10 GeV are rejected. 

The total visible invariant mass is required to be greater than 50 GeV/c 2 . 

The remaining events are then forced into two jets using the DURHAM-E algorithm. 
The polar angle of both jets, 6, is restricted to |cos6>| < 0.95. 



To reconstruct the effective centre-of-mass energy v s' , it is assumed that the ISR photon 
is emitted along the beam pipe, resulting in a boost of the // system, or the produced 
Z, in the opposite direction. A kinematic reconstruction based on rescaling of the jet 
energies is performed for each event. According to Monte Carlo studies, jet velocities and 
angles are well measured. Each jet four-momentum is scaled by a factor a according to 
the relations 
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— tree — wmis 

E\ ec = a x E~ ; P x = a x P x 
E r 2 ec = a 2 E~ ; P7 £C = a 2 P 2 meaS 



where rec stands for reconstructed and meas for measured. Correction coefficients are 
different for the two jets. 

In the absence of final state radiation (FSR), the correction coefficients a% and a 2 are 
obtained by applying energy and momentum conservation, and the assumption that one 
zero mass particle escapes detection (the ISR 7) along the beam axis. The solution kept is 
the one which yields ot\ and a 2 positive and closest to 1. Then, the effective centre-of-mass 
energy can be expressed as 

(pmeas pmeas \ 
tl ±2 
j^raeas ' ^meas I 

from which the quantity x = 1 — (s'/s) is reconstructed independently of the nominal 
centre-of-mass energy. 

The fit to the centre-of-mass energy is performed for events with 0.60 < x < 0.88 where 
the highest purity and the maximum sensitivity is achieved. The generated cross-sections, 
the effective ones after selection cuts and within the fitting window are summarised for 
each process in Table [3]. 



Processes 


tr(pb) 


o- e //(pb) 


^window (ph) 




108.84 ± 0.08 


68.61 ± 0.20 


40.07 ± 0.15 


WW 


16.02 ± 0.01 


10.51 ± 0.05 


0.47 ± 0.01 


zz 


2.545 ± 0.13 


1.08 ± 0.05 


0.33 ± 0.02 


77 — > hadrons (untagged) 


7800. ± 390. 


0.98 ± 0.22 


0.24 ± 0.11 


77 — > u/d 


474.00 ± 23.70 


0.61 ± 0.08 


0.19 ± 0.04 


77 — > ss 


26.00 ± 1.30 


0.04 ± 0.02 


0.01 ± 0.01 


77 — > cc 


93.60 ± 4.68 


0.55 ± 0.008 


0.19 ± 0.05 


77 — > bb 


0.53 ± 0.03 


0.02 ± 0.005 


0.01 ± 0.01 


Zee 


6.80 ± 0.27 


4.57 ± 0.18 


1.02 ± 0.04 


Wev 


0.608 ± 0.03 


0.32 ± 0.02 


0.11 ± 0.01 


Purity (%) 


78.60 ± 0.01 


93.93 ± 0.01 



Table 1: Generated (second column) and effective cross-section after selection cuts (third column) 
and within the fitting window (fourth column), for signal and background processes (first column) at 
V^=182.675. 



5 LEP Centre-of-Mass Energy measurement 

The LEP centre-of-mass energy is determined for each of the 181, 182, 183 and 184 GeV 
datasets. The four values are combined taking into account the relative luminosity of each 
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dataset. A Monte Carlo reweighting procedure [|14| is applied to find the value of Ecm 
which best fits the reconstructed x distribution. Selected Monte Carlo (KORALZ) signal 
events from the large sample at a reference energy of 182.675 GeV are reweighted using 



the ratio of the differential production cross-sections |TT5 



w ^ l cm) - 1 — i t^Jy 

~fa\ x truei ^CM) 

where Xi denotes the x at which the ith event has been generated. Background Monte 
Carlo samples are included in the fit, but they are not reweighted, i.e., their energy 
dependence is not taken into account. This introduces a systematic error which is 
estimated in section |7|. 

A maximum likelihood fit is performed with fixed bins of 0.01 over the x range of 0.60- 
0.88. The statistical error on the centre-of-mass energy is derived from the individual fits 
to the data distributions. 

The linearity of the reweighting technique was verified using five independent Monte 
Carlo samples (signal and background) generated within 2 GeV around the central 
nominal energy. The fitted values from these distributions have a slope 1.035±0.034, 
consistent with unity, and an insignificant offset of 0.045±0.049 GeV. 



6 Results 

The number of selected events at each nominal energy, and the corresponding number of 
expected events are summarized in Table ||]. The uncertainty on the number of expected 
events comes from three sources: the error on the integrated luminosity measurement, 
the statistical error on the estimation of the efficiencies and the theoretical uncertainty 
on the cross-sections. 



Nominal 

dataset 


Data 
events 


Expected 
events 


Integrated 
Luminosity (pb _1 ) 


181 


11 


8 ± 1 


0.166 ± 0.006 


182 


187 


173 ± 2 


3.924 ± 0.029 


183 


2293 


2191 ± 12 


50.795 ± 0.105 


184 


76 


84 ± 1 


1.927 ± 0.021 



Table 2: Number of events found in each energy dataset after the selection cuts, as well as the 
corresponding number of expected events and integrated luminosity . 



The LEP centre-of-mass energy values and the statistical errors obtained from fits 
to the individual datasets are summarized in Table |3], as well as the weighted average 
according to the integrated luminosity of each dataset. The last column gives the value 
obtained by the LEP energy calibration group. 

Figure [I] shows the reconstructed x distribution for the selected events in each dataset, 
superimposed with the Monte Carlo reweighted expectation using the LEP centre-of-mass 
energy which best fits the data. 
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Nominal 

energy 


EcM(flt) 

(GeV) 


-AEcm 
(GeV) 


+AE C m 
(GeV) 


E CAf (LEP) 
(GeV) 


181 


181.19 


—2.45 


+2.82 


"i on oo/" 1 i r\ r\ r - r\ 

180.826 ± 0.050 


182 


181.73 


-0.63 


+0.68 


181.708 ± 0.050 


183 


182.56 


-0.22 


+0.19 


182.691 ± 0.050 


184 


182.56 


-0.99 


+ 1.11 


183.801 ± 0.050 


Combined 


182.50 


-0.20 


+0.18 


182.652 ± 0.050 



Table 3: Fitted LEP centre-of-mass energies and corresponding statistical errors for each dataset. The 
last column gives the LEP centre-of-mass energies established by the LEP energy calibration group for 
the corresponding dataset. 
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Figure 1: Reconstructed x distribution for selected data (dots) and background (shaded area) for each 
energy dataset: 181 (a), 182 (b), 183 (c) and 184 GeV (d). The open histogram shows the results of the 
fit 
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7 Systematic uncertainties 



Table 0] summarizes all the systematic errors discussed below. 



kJW LLl CC 


AR„„ CMpV") 

alJ CM ^ivie v ) 


Calorimeter calibrations 


44 


Jet energy calibration 


47 


Jet angle calibration 


20 


MC statistics 


24 


ISR 


21 


ISR-FSR interference 


10 


Background contamination 


32 


Total 


82 



Table 4: Summary of the systematic errors on the LEP centre-of-mass energy measurement. 



Calorimeter calibrations 

During the 1997 data taking the uncertainties on the absolute calibration scale of 
ECAL and HCAL energies were estimated to be ±0.9% and ±2%, respectively. Maximum 
shifts of 26 MeV and 36 MeV were observed when applying the fitting procedure to 
the MC sample after variation of ECAL and HCAL energies in both directions by the 
corresponding amount. 

Jet energy calibration 

Z peak data are used to map the response of the detector to hadronic jets as a function 
of their polar angle. The observed difference in the energy scale of jets between data 
and Monte Carlo is parametrised as a function of |cos0|, where 9 is the angle between 
the jet direction and the beam axis. To estimate the systematic error, two modified 
parametrisations are derived which correspond to ±ler uncertainty on the discrepancy. 
The largest shift observed when these modified parametrisations are used to correct the 
jet energies amounts to 47 MeV and is taken as the systematic error. 

Jet angle calibration 

Z peak data have been used as well to study the two jets acollinearity as a function 
of |cos#|. The observed difference between data and Monte Carlo is used to estimate this 
systematic error. A 20 MeV shift is observed when applying or not this correction and is 
quoted as the systematic error. 

Initial State Radiation 

ISR is described in KORALZ up to second order in the leading-log approximation 
(0(a 2 L 2 )), in the YFS style. The effect of the missing higher order terms is studied 
by degrading KORALZ to 0(a 1 L 1 ). A 21 MeV shift is observed and is quoted as an 
upper limit of the systematic error. 

ISR-FSR interference 

No existing Monte Carlo generator is able to give a good representation of initial-final 
state interference for quark final states. An evaluation of the order of magnitude of the 
effect has been performed using ZFITTER []17||. The differential cross-section da / dx has 
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been computed with ZFITTER, switching on and off the ISR-FSR interference, including 
only the polar angle region |cos(0)| <0.95, at the energy corresponding to the fitted value 
in the data. In both cases these differential cross-sections are integrated in x bins of 0.01 
over the x range of 0.60-0.88, as for the energy measurement. The shift observed in the 
Z resonance peak translates into an uncertainty of 10 MeV in the determination of the 
LEP centre-of-mass energy. 

Background contamination 

The expected background level after selection cuts is 6% and the uncertainty on 
the contribution of the different subprocesses will affect its shape which may also 
depend on the centre-of-mass energy. Since the background is not reweighted, this 
introduces an additional source of systematic error both from the background shape and 
its normalization. 

Background shape. The normalization uncertainty is estimated independently for 
each background process. 



The two-photon background, simulated with PHOT02 and PHOJET, is normalized 
to the data in the low visible mass region, M vis < 50 GeV/c 2 ; the difference between 
the expected and normalized cross-sections is used to estimate the uncertainty of 
the process, as in [|Kj . 



For the WW, ZZ and Zee processes the uncertainties are computed from the relative 
variation of their respective total production cross-sections from the generated 
energy (183 GeV) down to the fitted energy from the data. This is calculated using 
GENTLE m for the WW process and PYTHIA for the ZZ and Zee process. 



Each background process is varied according to its uncertainty keeping constant the 
overall background level. The largest shift observed in each case is quoted in Table |^, 
leading to a total systematic error of 8 MeV. 



Process 


Deviation (MeV) 


WW 


4 


ZZ 


2 


77 


6 


Zee 


3 


Total 
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Table 5: Systematic error on the centre-of-mass energy due to background shape uncertainties. 



Total normalization. The systematic error coming from the uncertainty in the total 
background normalization is estimated by varying all the background processes up and 
down, according to their respective uncertainties. The largest shift observed is 31 MeV 
and is quoted as the systematic error. 
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8 Summary and conclusions 



The LEP centre-of-mass energy can be determined from the kinematic reconstruction of 
qcpi events with a qq invariant mass around the Z mass. The average LEP centre-of-mass 
energy at ALEPH for the high energy run of 1997 is measured to be 

E CM = 182.50 ± 0.19 (stat.) ± 0.08 (syst) GeV 

This result is consistent with the estimate from the LEP energy working group f|| 
E CM = 182.652 ± 0.050GeV. 

With the expected increase in statistics at LEP2 and with refined experimental 
techniques akin to those used for the W mass determination the method described 
here should provide an alternative measurement of the LEP centre-of-mass energy, with 
competitive accuracy. Already, if only this evaluation of the LEP centre-of-mass-energy 
were used in the W mass determination, the systematic error on Mw coming from the 
precision on the beam energy scale would be 90 MeV, (AM W / M w = AEcm/Eqm), which 
is smaller than the experimental error of 139 MeV on My/ from the same data sample 0. 
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